Magnetic field dependence of spin-lattice relaxation in the s± state of Bao.67Ko.33Fe2As2 
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The spatially averaged density of states, {N(0)), of an unconventional d-wave superconductor 
is magnetic field dependent, proportional to H^^^, owing to the Doppler shift of quasiparticle ex- 
citations in a background of vortex supercurrents.— This phenomenon, called the Volovik effect, 
has been predicted to exist for a sign changing s± state;^ although it is absent in a single band 
s-wave superconductor. Consequently, we expect there to be Doppler contributions to the NMR 
spin-lattice relaxation rate, l/Ti oc{A*'(0)^), for an s± state which will depend on magnetic field. 
We have measured the ^^As l/Ti in a high-quality, single crystal of Bao.67Ko.33Fe2As2 over a wide 
range of field up to 28 T. Our spatially resolved measurements show that indeed there are Doppler 
contributions to l/Ti which increase closer to the vortex core, with a spatial average proportional 
to , inconsistent with recent theory— 

PACS numbers: 



The spin-lattice relaxation rate, 1 /Ti , in the supercon- 
ducting state can provide valuable information about gap 
structure and about the effects of vortices on the quasi- 
particle density of states. There have been a number of 
reports on the temperature dependence of 1/Ti at low 
magnetic fields for various pnictide superconductors)^^— 
which have been interpreted in terms of order parameter 
structure. However, impurities of unknown origin and 
concentration can play an important role in the analysis. 
For this reason, conclusions about the superconducting 
state based on temperature dependences can be ambigu- 
ous. 

An alternative approach to explore unconventional 
character of the order parameter is the magnetic field de- 
pendence of the density of states, which can be specific 
to a particular order parameter symmetry, easily probed 
through specific heat or thermal conductivity measure- 
ments?^ The unconventional structures of the s± and 
d-wave states each lead to characteristic magnetic field 
dependences of the spatially averaged density of states, 
attributable to the Doppler shift of quasiparticle exci- 
tations, known as the Volovik effect ^i— For the d-wave 
case c-axis line nodes in the gap give rise to a field de- 
pendence, cx H^/^ and a non-linear Mcissncr effcctJ^ 
In contrast, according to Bangp. in the case of s± sym- 
metry for a multiband superconductor, the spatially av- 
eraged density of states at the Fermi energy is propor- 
tional to the magnetic field. Since 1/Ti is proportional 
to the square of the local density of states, one might 
think that the predicted Volovik effect should be oc H^. 
However, the Volovik effect pertains to the spatial av- 
erage over the vortex unit cell which decreases in area 
inversely proportional to the magnetic field. According 
to the theory, there is a region of normal state excitations 
surrounding the vortex core of radius ^(A2/Ai), where 
£, is the core radius (~ 30 A), equal to the coherence 
length, and A2/A1 is the ratio of large to small gaps, 
leading to the prediction, \/TiT (x H at low magnetic 



fields, H < i7c2(Ai/A2)2 - ST; otherwise it should 
be constant^ In this Letter, we report ^^As NMR mea- 
surements in single crystals of Bao.67Ko.33Fe2As2 cover- 
ing a wide range of magnetic fields. Our results show 
that indeed there is a Doppler contribution to the spa- 
tially averaged spin-lattice relaxation but that in the low 
temperature limit, l/TiT oc over the whole range of 
magnetic field, inconsistent with predictioni^ 

We performed our ^^As NMR measurements at North- 
western University and the National High Magnetic Field 
Laboratory, from 4 K to room temperature with external 
magnetic field from 6.4 to 28 T. The fields were paral- 
lel to the c-axis of the single crystals, Bao.67Ko.33Fe2As2 
(BaK122) that had a zero-field Tc = 38 K and were grown 
at the University of Tennessee by the self-fiux method 
To increase signal intensity in the superconducting state, 
the crystals were cleaved to dimensions of 3x3x0.1 mm'^ 
and total mass of 17 mg. Typically, spin echo sequences 
{tt/2 - tt) were used to obtain the spectrum. Knight shift, 
and 1/Ti for the central transition (-1/2 4=^ 1/2) with a 
TT-pulse ~ 7 /isec. The spin-lattice relaxation was mea- 
sured with the full recovery method (28 to 300 K) and 
progressive saturation techniqueSfi^ (4 to 26 K) the latter 
being more accurate for very long relaxation times at low 
temperatures. The average rate was measured with the 
TT-pulse centered on the spectrum. Frequency-resolved 
spin-lattice relaxation was also measured by dividing the 
spectrum into many small frequency windows and the 
relaxation was determined separately in each window. 
Knight shift measurements were performed with a fre- 
quency sweep method. 

Early experiments on optimally, hole-doped, single 
crystals of Bao.6Ko.4Fe2 As^^ grown with tin ffux did not 
detect any signal below 20 K due to linewidth broaden- 
ing from paramagnetic impurities on the As sites at a 
level of 1%. However, there have been substantial 
improvements in lowering the impurity concentration us- 
ing the self-flux method.W.'i3. The frequency-swept ""^As 
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FIG. 1: As NMR spectra of Bao.67Ko.33Fe2As2 measured 
by a frequency sweep technique in 13 T with //He-axis of 
the crystals. Below Tc = 32 K (blue trace) the spectra shift 
noticeably to lower frequency with decreasing temperature. 
The linewidths of the spectra decrease in the superconducting 
state as reported in some other compounds,— where this 
was attributed to reduction in the local field distribution from 
impurities. 



NMR spectra of our crystals in 13 T with ff||c-axis, are 
shown in Fig. [H The in H = 13 T is 32 K, and 
a shift of the spectra can be easily seen. This decrease 
of the Knight shift indicates spin-singlet pairing in the 
superconducting state. On cooling, the linewidth slowly 
broadens from 60 kHz at T = 300 K to 70 kHz at Tc. 
Below Tc the linewidth increases up to 110 kHz near 20 
K, and then it decreases to 80 kHz at 4 K, and is inde- 
pendent of magnetic field from 6.4 to 16.5 T, to within 
10%. The weak dependence of the linewidth on mag- 
netic field and temperature in the normal state indicates 
that few magnetic impurities are present, comparable to 
the cleanest cuprate crystals such as Bi2SrCa2Cu2 08+i5 
(Bi2212)i^ This point is also consistent with the similar 
results we find from our comparison of the zero field ex- 
trapolations of 1/Ti with those of clean Bi2212 crystals 
which we discuss later. The Knight shift, K = Kg + Korb, 
was determined from the first moment of the NMR spec- 
trum where Kg and Kort sltc the spin and orbital parts 
of the shift, respectively. The orbital part is temperature 
and field independent, consequently the temperature de- 
pendence of the shift in Fig. [5] can be associated with 
Kg, decreasing below Tc on cooling. The solid curve in 
the figure is the temperature dependence of Kg that we 
describe with a phenomcnological model for the density 
of states, Eq. [1] based on the parameters obtained from 
1/Ti measurements. 

The behavior of the spin-lattice relaxation in the su- 
perconducting state is the main focus of our present work 
where we measure the temperature and magnetic field de- 
pendence for H = 6.4, 10.8, 14, 16.5, 27 and 28 T, parallel 
to the c-axis of the crystals. The rates were measured 
with the spectrometer frequency set at the peaks of the 
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FIG. 2: The total Knight shift K(T) is shown for H ^ 13 T 
where the temperature dependence is associated with the spin 
part, Ka, that decreases below Tc, consistent with spin-singlet 
superconductivity. The black arrow indicates Tc = 32 K. The 
data can be fit phenomenologically assuming that the low 
temperature spin shift is proportional to an average density 
of states, Eq.[T] represented by the solid curve, with an orbital 
shift, Korb = 0.21%. 



spectra. A coherence peak below Tc was not observed, 
and the suppression of Tc by the magnetic field was min- 
imal, from T = 32 to 30 K when the external field was 
increased from 6.4 to 27 T. In low magnetic fields, i.e. 6.4 
and 10.8 T, the temperature dependence of l/Ti could 
be approximately described as T^ at intermediate tem- 
perature, as has often been reported elsewhereAii But 
in higher fields, 14, 16.5, 27 T, below T = 10 K, we find 
1 /Ti cx T , indicating a constant average density of states 
at zero energy , (iV(0)). Recently Li et ali^ observed 
an exponential temperature dependence of the rate in a 
magnetic field oi H = 7.5 T consistent with the presence 
of a full gap. A comparison of our data with that of Li 
et al. shows that they are identical except at our lowest 
temperature, T = 4 K, where our higher value of l/Ti 
might be understood as the effect of residual impurities 
in our crystal obscuring exponential behavior. Increasing 
the magnetic field we find that the spin-lattice relaxation 
at 4 K increases systematically indicating the existence 
of a field dependent density of states at the Fermi sur- 
face. This observation is a characteristic signature of a 
Volovik effect. 

We use a phenomcnological model to fit 1/Ti in vari- 
ous magnetic fields. We express the thermal and spatial 
average over the density of states at the Fermi surface as, 

(7V(0)) = a{H) + 6oe-^^/'=«^ + cqc-^^/'^^^ (1) 

where a[H) ~ ao + a-iH + .. and gq represents possi- 
ble contributions from non-magnetic impurities. The two 
gaps, Al and A2, appear in exponential terms with rela- 
tive weights, bo and Cq, as might be expected for the low 
temperature limit. Since 1/TiT oc (A^(0)^), our model 
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FIG. 3: The spin lattice relaxation rate, 1/Ti, of 
Bao.67Ko.33Fe2As2 in magnetic fields, H = 6.4, 10.8, 14, 16.5, 
and 27 T. The rates were measured at the peak position in 
the spectrum. The data is consistent with our two-gap model 
(solid curves) provided the magnetic field dependence at low 
temperature is oc . The black arrow indicates Tc- 



for 1/TiT becomes, 

1/TiT cx [a{H) + bae-^''^^^ + coe-^'''''^'^f. (2) 

At low temperatures the two exponential terms arc of lit- 
tle importance and the rate is determined by a{H). Our 
numerical analysis provides fits for all of the parameters 
of the model. Below H = 16.5 T, we take them to be 
magnetic field independent. However, at this and higher 
magnetic fields we find that the relative weight of the ex- 
ponential term from the smaller gap, bo, must be reduced 
compared to the larger gap weight, cq, in order to fairly 
represent the data. As stated previously, these gap pa- 
rameters are not important in the low temperature limit 
where we seek to describe the field dependence of the re- 
laxation rate and so we do not ascribe specific importance 
to this additional field dependence other than it allows us 
to represent the high temperature behavior in each field. 
Nonetheless, we point out that our results for the tem- 
perature dependence at low magnetic field are identical 
to those from Li et ali^ for clean crystals, except for the 
lowest temperature point at 4 K. 

Our analysis in each field is shown in Fig. IHa)-(e), 
where Ai and A2 are 2.1 ± 0.2 meV and 12.1 ± 1.4 meV 
respectively. The sizes of the gaps correspond well to 
the sizes of the 3D superconducting gap function from 
ARPES measurements, 2.07 meV and 12.3 meV.^-S The 
ratio of the coefficients, b^ and cq, decreases ai H = 16.5 
and 27 T, indicating a possible suppression of the smaller 
superconducting gap, Ai, by the external magnetic field. 
The low temperature magnetic field dependence of 1/TiT 
is given by a(H) shown in Fig. HJ^f). The behav- 
ior might be associated with Doppler shifted quasiparti- 
cles, although the field dependence is different from that 
predicted by theory^ It should be noted that the elec- 
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FIG. 4: (a)-(e) l/^^TiT with the best fit curve in each mag- 
netic field. The unit of the a(H), bo, and co is (s~^ K~^), and 
the gaps are in meV. (f) l/^^TiT in BaK122 at 4 K, is shown 
for comparison with the Zeeman contributions to l/^^TiT in 
Y123 and Bi2212 at 5 K. Assuming that the electronic (j- factor 
is the same for BaK122 as for the cuprates, we argue that the 
Zeeman contribution to the field dependence of the average 
rate we have measured in BaK122 is significantly smaller than 
from Doppler contributions. In the case of the cuprates, the 
Zeeman contributions to the spin-lattice relaxation rate were 
isolated using frequency resolved measurements performed at 
the saddle point of the local field distribution (peak of the 
spectrum) where Doppler contributions cancel based on sym- 
metry of the supercurrents from near-neighbor vortices. 



tronic Zeeman interaction also contributes to the quasi- 
particle energy giving a H"^ dependence to 1/TiT^ In 
Fig. mi) we show the field dependence of l/TiT for ^'^O 
NMR from YBa2Cu307+5(Y123) ahgned powdersiS and 
Bi2212 crystals^ which has been attributed to this Zee- 
man term. From comparison with these compounds, al- 
lowing for the 27% larger gyromagnetic ratio of arsenic 
compared to oxygen, it is reasonable to conclude that 
the significantly larger field dependence of \/TiT for 
BaK122, cannot be attributed to the Zeeman term. It 
is notable that the three materials have a similar value 
of the spin-lattice relaxation in the limit of zero field. 
Although, to some extent, this might be fortuitous, it 
nonetheless suggests that our BaK122 crystal does not 
have significantly more impurity scattering than these 
high quality cuprate materials. Since the NMR linewidth 
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at r = 40 K is independent of magnetic field for H < 16.5 
T to within 10%, we do not associate the field depen- 
dence of the rate with magnetic impurities. However, 
this possibility can be investigated further by measure- 
ment of the frequency-resolved spin-lattice rate which we 
describe next. 

For unconventional superconductors 1/Ti can depend 
on the position of the probe nucleus relative to the vor- 
tex coreii^Ti^i The increase in the supercurrent momen- 
tum, Ps, approaching the core leads to a corresponding 
increase in the Doppler shift of the energy of quasiparti- 
cle excitations, vp ■ Ps, where vp is Fermi velocity. The 
vortex core, having the highest local magnetic field, cor- 
responds to the largest frequency in the NMR spectrum. 
We have looked for evidence of this spatial dependence of 
l/Ti through frequency-resolved, i.e. spatially resolved, 
measurements performed across the spectrum, as shown 
in Fig. [SJ 

In the normal state (40 K) we find a fiat l/Ti distribu- 
tion throughout the spectrum as expected in the absence 
of Doppler terms or magnetic impurity contributions to 
the rate. In the superconducting state, there is an in- 
crease of 1/Ti with frequency, developing markedly at 
T = 26 K with more than an order of magnitude varia- 
tion across the spectrum. 



even in a somewhat disordered vortex structure, the high 
field portion of the spectrum can be associated with nu- 
clei in the vortex core. This is the case for the distribu- 
tion in 1/Ti observed in Y123, which was attributed to 
the Doppler shifti^ of quasiparticle energy from vortex 
supercurrents. Our frequency-resolved measurements of 
l/Ti in BaK122, Fig. 5, show the existence of a spatially 
inhomogeneous distribution which onsets with supercon- 
ductivity. We ascribe this to the vortex state for which 
the most likely explanation is a Volovik effect. Another 
explanation was suggested some years ago to explain ob- 
servations in superconducting vanadium compoundsi^i^ 
There it was argued that spin-diffusion from relaxation 
sources in the vortex core might produce a spatially in- 
homogeneous distribution of l/Ti. Later measurements 
and theoretical work by Genack and Redfiel d^^i^^ showed 
that this suggestion was incorrect, and that spin diffusion 
is quenched on very short time scales owing to depiction 
of the dipole energy reservoir, an effect even further sup- 
pressed with increasing field. We measured the spin lat- 
tice relaxation rates in higher fields, 24 T and 28 T, as 
shown in Fig. [B) An inhomogeneous spin-lattice relax- 
ation rate distribution was found similar to that of H = 
16.5 T, Fig. [5l and rules out spin diffusion as a possible 
mechanism ,2i2i 
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FIG. 5: Spin-lattice relaxation rate across spectrum in the 
normal state (a) and superconducting states(b),(c),(d) in 16.5 
T, with _ff||c-axis. In the normal state, at 40 K, there is no 
significant frequency dependence in l/Ti. However, the rate 
becomes dependent on frequency as the sample is cooled deep 
into the superconducting state. 

We note that the linewidth, ~ 80 kHz at 4 K in 
H = 16.5 T is somewhat broader than our calcula- 
tion from Ginzburg-Landau theory using Brandt's algo- 
rithm^ for a perfect vortex lattice, ^ 23 kHz. However, 




0.01 
174.9 175.0 




0.01 



175.1 175.2 175.3 204.1 2042 2043 204.4 204.5 
Freq.(MHz) Freq.(MHz) 



FIG. 6; Spin-lattice relaxation rate across spectrum at 4 K 
in 24 T (a) and 28 T (b) with _ff||c-axis. The inhomogeneous 
frequency dependence of the rate is observed in both mag- 
netic fields similar to H = 16.5 T, Fig[5] (d). Additionally, 
the increase in the high frequency part in the spectra can be 
understood as an asymmetry from the vortex lattice. 



With reports from experiments in cuprates a decade 
g^gQ i9,27 ^Yiis mechanism was studied theoretically by 
WortiS)^ who came to the same conclusion. A more 
detailed discussion has been provided by Mounce et al^ 
We point out that in the recent theory^ii of the Volovik 
effect in s± superconductors the combined effects of the 
Zeeman interaction and vortex supercurrents have not 
been taken into account. Their importance was indicated 
in the work of Mitrovic et al. on YBa2Cu307^ and might 
be an important component missing from the theory. We 
conclude that our observations are most likely a conse- 
quence of vortex supercurrents but for which there is not 
yet a satisfactory theoretical explanation. 

In summary, we have studied the ^^As Knight shift 
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and spin-lattice relaxation rate in slightly underdoped 
Bao.67Ko.33Fe2As2 crystals in the superconducting mixed 
state. We found that 1/TiT approaches a constant at low 
temperatures in high magnetic field and is proportional 
to the square of the field. Although this is inconsistent 
with a theory for the Volovik effect,— our results can 
be accounted for by a phenomenological model which is 
based on s± symmetry with two isotropic gaps, and non- 
magnetic impurities. The distribution of 1/Ti across the 
spectrum resembles that observed in a vortex solid of an 
unconventional superconductor associated with spatially 
resolved Doppler contributions to the quasiparticle exci- 



tation spectrum. 
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